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Abstract 

We considered the dynamics of gravitating gas - a continuous media 
with peculiar properties. The exact solutions of its Euler equations 
for simple initial conditions is obtained. 

1 Introduction 

Gravitating gas is a mechanical system which could be considered as a con- 
tinuous limit of the iV-particles system described by Hamiltonian: 

R _Y^Pi 2 j2 GniiTn i ^) 

i=l 2 m i i^j \ x i ~ x j\ 

where Xi,pi are canonical coordinates of particles (stars) with masses rrii, 
G-gravitational constant. 

When N — > oo(with assumption rrii = m) the system turns into the 
kind of continuous media with very peculiar properties due to attractive 
interaction between constituents. In the usual gas the effective interaction 
between particles (atoms or molecules) is repulsive and as a result being put 
in a volume it spreads in space, filling after some time the whole volume with 
uniform density. On the opposite, as was pointed out in [T], the gravitating 
gas, because of attractive interaction may form isolated steady configuration 
with asymptotically vanishing density. In particularly, the galaxy named 
Hoag's object may be an example of such gravitational soliton [lj. Also, 
as we shall see below there is no sound waves in gravitating gas and local 
perturbation of density in linear approximation creates instability. That is 
why in this system in order to study evolution of perturbation we need to 
take into account nonlinearity of the equations of motion. In the present 
paper we shall consider the cases, where we can solve equations of motion 
exactly. 
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2 Equations of motion 

In fluid (gas) mechanics there are two different pictures of description. The 
first, usually refereed as Eulerian, uses as the coordinates the space dependent 
fields of velocity and density. The second, Lagrangian description, uses the 
coordinates of the particles x(£i,t) labeled by the set of the parameters £j 
(the numbers of particles in (pQ)), which could be considered as the initial 
positions £ = x(£i,t = 0) and time t. The useful physical assumption is that 
the functions t) define a diffeomorphism of R 3 and the inverse functions 
£(Xi,t) should also exist. 



The density of the particles in space at time t is 

p(x, t)= J d^p {^)5{x - x i)), (3) 

where po(0 is the initial density at time t — 0. The velocity field u as a 
function of coordinates x and t is: 

u(xj,i) = x(f(x i9 t),*), (4) 

where £(x,t) is the inverse function (j2J). The velocity also could be written 
in the following form: 

,s = J d 3 £po(£ i )x(£ i ,t)5(x- x(£j,t)) 

or 



p(xi,t)v(xi,t) = j d£po(tii)x(€i,t)8(x-x(€i,t)). (6) 
Let us calculate the time derivative of the density using its definition (J4]) : 

p(xi,t) = J d 3 £p (£i)—5(x- x(£i,t)) 

= J d 3 Mtii) (-£(&,*)) 

= -i / ^ 3 epo(6)^(6, - x(&, t)) 

d 

= -—p(xi,t)v(x h t) (7) 



In such a way we verify the continuity equation of fluid dynamics: 

p(xi, t) + d(p(xi, t)v(xi, tfj = 0. (8) 



Using the coordinates x(£i,t) as a configurational space variables we can 
write the Lagrangian for the continuous generalization of the system, de- 
scribed by (pQ) 



L = / d 3 £po(& 



mx t) Gm 



+ 



d 3 £d 3 £' 



Po(&)po(£0 



\x(&,t) - x(&t)\' 



(9) 



The equations of motion, which follow from the Lagrangian ( TTUj) have the 
form: 



mx(£i,t) + Gm / d £ Po(£i 



0. 



(10) 



Now we need to translate the equations (flUI) into the language of Euler 
variables. For this let us differentiate both sides of the equation ([6]) with 
respect to time: 



d t 

—p(xi,t)v(xi,t) = J d 3 £p {£i)x(£i,t)8(x - x(£i,t)) 

d 



(11) 



Substituting x(£i,t) from the equation ffTUj) and transforming the second 
term, as we did in derivation of the equation (j7])we obtain: 



^p(x h t)v(xi,t) 







(p(Xi,t)v(Xi,t)v k (Xi,t)) + 



^Po(6) -7 / d 3 C P o(i 



dx k 



5(x-x(£i,t)) 



(12) 



where we introduced notation 7 = Gm. To transform the last integral we 
first perform integration over the £: 



tfWfc) -7 / rf 3 e'po(c: 



7 p(^,t) / d 3 eVo(e; 



5(f - £(&,*)) 



|f-f(4',t)r 



(13) 



Now let us insert in the integral over £' the unity 



1 = J d 3 y8(y-x(&t)) 



(14) 



and change the order of integration: 




(15) 



where in the last step we have used the definition of p(yi,t) given by (j3J). 
Finally the equation of motion (ITT]) in terms of Euler variables takes the 
following form: 



Using the continuity equation we can rewrite (ITS"]) in the more familiar form 
of Euler equation: 



3 Absence of sound waves in gravitating gas. 

As very well known (see for example [2]), small perturbation of density in gas 
or fluid (and/or velocity) creates the sound waves. In the case of gravitating 
gas, due to attractive interaction of constituents the situation is completely 
different. Let us introduce a small perturbation of density p(xi,t), so that 



where po(xi) is background density. Neglecting nonlinear term in Euler equa- 
tion we have 




(16) 




(17) 



p(xi,t) = p + p(xi,t), 



(18) 




(19) 
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In linear approximation the continuity equation will have the following form: 

p{xi,t)+po — — — = 0. (20) 

Taking time derivative of equation f|20|) and substituting v(xi,t) from equa- 
tion f|T9|) we obtain: 

ftxi,t) + mA [d 3 y^^r = 0. (21) 
J \x — y\ 

Action of Laplace operator on the potential gives the density times — 4tt and 
finally we arrive at 

p{x h t) - ATT Pol (p (xi) + p{x h t)) = 0. (22) 

The solution of this equation has exponential time dependence which makes 
an assumption of small perturbation invalid. Also this example shows that 
the linear approximation does not work for gravitating gas. The reason 
for such behavior of perturbation is the attraction forces acting between 
constituents. In the case of usual gas with repulsive interaction between 
constituents the r.h.s. of the equation ffT9l) will stay e.g. ~~~ 8p ^g'^ , where 
k is the velocity of sound and therefore we will obtain the instead of (22) the 
usual wave equation: 

p(x i ,t) + k?Ap(x i ,t))=0. (23) 

The conclusion from this observation is that if we want to study the behavior 
of gravitating gas, the linear approximation is not sufficient and we need to 
solve Euler equations exactly. 

4 Exact evolution of spherically symmetric 
configuration of gravitating gas 

Here we will consider the spherically symmetric configuration of gravitating 
gas which is described by p(xi,t) = p(r,t) and v(xi,t) = -v(r,t). Substitut- 
ing this parametrization into equations (IE]) and ( JT71) we obtain: 

v(r, t) + v(r, t)d r v(r, t) = 7 <9 r / d*x'^§-, (24) 
r 2 p(r, t) + d r (r 2 p(r, t)v(r, t)) = 
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As is well known for spherically symmetric mass distribution the force in the 
r.h.s. of the first equation (25) is defined by the total mass inside the sphere 
of radius r: 

7*/dV^ = -W?M, (25) 
J \x — ar I r z 

where we introduced notation 

m(r,t) = f dr'r' 2 p(r\t) (26) 
Jo 

With this notation we can present the equations (24) in the following form: 

v(r, t) + v(r, t)d r v(r, t) = -4^^ (27) 
rh(r, t) + v(r, t)d r m(r, t) = 0. 

Let us consider the function f(r, t) which is defined by the following equation: 

f(r, t) + v(r, t)d r f(r, t) = 0, f(r, 0) = r. (28) 

The function m(r, t), because it satisfies the same equation could be written 
in the following form: 

m(r,t)=m(f(r,t)), (29) 

where m(r) is the initial data for function m(r,t): m(r) = m(r, 0). In such 
a way we can rewrite the first of equations (27): 

v(r, t) + v(r, t)d r v(r, t) = -Aix ^^^ ^ . (30) 

Multiplying both sides of this equation by v(r,t) and taking into account 
equation (28) we arrive at the following statement: 

(a + ^r,*)ft)[^-4» 7 =!M] = 0. (31) 

From this statement follows that the expression in the brackets is a function 
of f(r, t) only. So, we have 

-J^l_ 4nJ ^fMl = AWrA) . (32 ) 
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Putting in (32) t = we can express the unknown function A(r) via initial 
data for v(r,t) and m(r,t): 

v(r, 0) 2 m(r) , , . , . 

-±^-4n 1 ^ = A(r). (33) 

Therefore we obtain 

v(r,t) 2 v(f(r,t)) 2 A ... ,, r l 1 . 

where we put v(r, 0) = v(r). Apparently from (34) we can find the velocity: 
v(r,t) = [2 C U{r : t)f + 4^(/(r,t))[- - t^])] 172 . (35) 

This expression for velocity has to be substituted into equation (28) for func- 
tion f(r, t) which we need to solve. In order to do it let us introduce the 
function s(r, t): 

r = f(r,t)-s(f(r,t),t). (36) 
It is obvious that s(r, 0) = 0. Also from this definition follows that 

v{r,t) = -dts{f,t). (37) 

where we deliberately suppressed the arguments of the function f(r,t) to 
make it clear that the differentiation in (37) is with respect to the second 
argument of s(f,t). (This relation could be obtained by applying operator 
d t + v(r, t)d r to both sides of (36)). In such a way we have 

*,( f\2 i 1 

dts(f,t) = -[2(-^- + 4 7 r 7 m(/)[ / _ g(/ ^ ) - -])] 1 / 2 . (38) 

This relation gives us a possibility to find time dependence of the function 
s(f,t): 

r-(/.*) ds _ 

o [2(^f+4 7 r 7 m(/)[ 7 i i -}])]V2 

Knowing the function s(f,t) we can find function / making use of Lagrange 
formula for analytic branch of inverse function In the present context 
this formula takes the following form (here we will suppress the argument t of 
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the function s(f,t), considering it as a parameter). So we have the following 
relation: 

r = f - s(f) (40) 
The inverse function f(r) is given by 

f = r + s(r) + l -{s\r))' + ^(s 3 (r))" + ^(s\r))^ + ■■■ + ^(s»)^ 

(41) 

This formula gives the complete solution of our problem, though in that form 
it is not at all clear what is the evolution even for the simplest cases. As 
an example we shall consider the case where we can at least qualitatively 
describe evolution. The first simplification we shall do is to put v(f) = 0. In 
this case the equation (39) takes the following form: 

- { S, lmU )r^ C ^r^-t. (42, 

Now we shall change the integration variable (this parametrization is usual 
for Kepler problem [I]) 

, = /(1 (43) 



after which we obtain 



$7T7m(/) 



^ L a + sina] = t. (44) 



Next simplification consists in choice of the function m(r). If we take uni- 
form initial density in the whole space p(r) = Pq, then m[r) = Poy and 
equation(44) becomes 

3 I 1/2 1 . 
~\a + sina\=t. (45) 

Let us summarize the result. The equation (36) gives us the relation 

r = /i±p. (46) 

Time dependence of a is given by equation (45), which we can rewrite in the 
following form 

t . . 

a + sina = —2ir—, (47) 



where we introduced notation T = ty 



1/2 



87T7P0 . 



Note that in this case we 



don't need to use Lagrange formula in order to find the function f(r,t). 
From equation (46) we have 



f(r,t) 



2r 



1 + cosa{t) ' 

from where we obtain the evolution of the functions m(r, t) and v(r, t): 
m(r, t) = m(f(r, t)) 



(48) 



2r 



_ l+cosa(t) 

„,( r f\ _ d t f _ r 87r7p il/2 2rsina(t) 
u \'i l ) d r f L 3 J (l +COSQ (t))' 

Also from m(r, t) we obtain the evolution of density: 
p(r,t) = \d r m(r,t) = p 



1 + cosa(t) 



(49) 



(50) 



In all these equations time dependence of a(t) is defined by equation (47). 
Note that according to it a(t) is zero at t — and moves in negative direction. 
We see the the density is unform for all t, oscillating between p and infinity. 



5 Evolution of flat configuration gravitating 
gas. 

Let us consider the case where v = (v(x, t), 0, 0) and p = p(x,t). Euler 
equations takes the following form: 

v(x, t) + v(x, t)d x v(x, t) = —2ix^d x J dy\x — x'\p(x', t) 

p(x,t)+d x (p{x,t)v(x,t)) = (51) 

This form could be obtained from the 3D equations and the kernel — 2ix\x— x'\ 
is the result of integrating rgzgr, over variables which don't enter into p. 
Needless to say that this one-dimensional gravitating gas could be interpreted 
as flat layer in three dimension. Differentiating potential in the r.h.s of Euler 
equation we can rewrite it in the following form: 

v(x, t) + v(x, t)d x v(x, t) = — 2tt7 / dye(x — x')p(x', t) 

= -2^7^ dx'p(x', t) - r dx>p(x', t)}, (52) 
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where e(x) = sign(x). It will be convenient to introduce instead of density 
another variable g(x,t) 

If 1 r x f°° 

g(x,t) = — I dye(x — x)p(x ,t) — -(/ dx'p(x',t)— / dx'p(x',t)). (53) 

Z J Z J— oo J x 

Apparently this new variable plays the same role as the function m(r, t) in 
the previous section. With this new variable the the system ( IBTT) will take 
the following form: 

v(x, t) + v(x, t)d x v(x, t) = —47T'jg(x, t) 

g(x,t) + v(x,t)d x g(x,t) = 0. (54) 

Now let us introduce the first integral of the second equation f(x,t), satis- 
fying initial condition f(x, 0) = x. Apparently, the function g(x, t) could be 
written as 

g(x,t) = g(f(x,t),0), (55) 

so, evolution of g(x,t) is just reparametrization of its initial data g(x,0) = 
g(x). Now it is easy to find general solution of the first equation (1541) : 

v{x, t) = v{f{x, t)) - t4n ig (f{x, t)), (56) 

where v(x) is the initial data for velocity. The first term in f)56p is the solution 
of homogenous equation which we fixed through the initial data for v(r,t), 
while the second is the solution of inhomogeneous one. Now we have arrived 
at the point similar to that in the previous section. We need to find the 
function f(x,t). We shall do it as before. Let us introduce the function 
s(f, t) through the equation 

x = f(x,t)-s(f(x,t),t), (57) 

such that 

d t s(f,t) = -v(x,t), (58) 

where the derivative is taken with respect to the second argument. The case, 
we are considering now is simpler when the previous one and we can find 
s(f, t) explicitly: 

s(f,t) = -tv(f)+^ t -g(f). (59) 
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Now we again can use Lagrange formula to find f(x, t) 



f = x + s(x) + l -{s\x))' + ^(s 3 (x))" + ^(s\x))^ + ■■■ + i(^)) (B_1) 

(60) 

Remember that here we suppressed the second argument of the function 
s(f,t), considering t as a parameter. 

So, formally the Euler equations are solved if we can evaluate for given 
initial data p (x),v(x) the function f(x,t) presented by infinite series fl60|) . 
Again, for some simple initial data it is possible to solve equation fl59|) without 
Lagrange formula. For example let us take v(x) = and density po(x) = b 
for \x\ < a and po(x) = outside this interval. In this case the function s(f) 
in (15 9p is linear and we can find explicitly evolution of gravitating gas. 

{—ab, x < —a(l - 27vjt 2 b), 

_^_ ) -a(l-2n 1 t%) <x<a(l-2^ 7 f 2 6), (61) 
ab, x > a(l - 2-K^t 2 b). 

Differentiating g(f(x,t)) by x we obtain from (1M|) evolution of density 

TO, x < -o(l -2^ 2 6), 

P(^> *) = | (i-2^6) » - a ( X " 2yr ^ 2 &) < x < a(l - 27r 7 t 2 fo), (62) 
I 0, x > o(l - 27r 7 t 2 6) 

and from (|58|) we can get also evolution of velocity. This solution of Euler 
equations which we o btained describes collapse of gravitating gas during 
finite time T = to zero size. Here arises the question what happens 

after this time. The analysis of the Lagrange equation shows that the motion 
continues in backward evolution. Gas emerges from one point and returns 
to the initial configuration. The impossibility of analytic description of the 
whole evolution has the same origin as in the problem of bouncing ball. 
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